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ABSTRACT 

We study in simple terms the role of feedback in establishing the scaling relations 
of low-surface-brightness (LSB) and dwarf galaxies with stellar masses in the range 
6 x 10 5 s$ M* ^ 3 x 10 10 M Q . These galaxies, as measured for example from SDSS 
and in the Local Group, show tight correlations of internal velocity, metallicity and 
surface brightness (or radius) with M». They define a fundamental line which distin- 
guishes them from the brighter galaxies of high surface brightness and metallicity. The 
idealized model assumes spherical collapse of CDM haloes to virial equilibrium and 
angular-momentum conservation. The relations for bright galaxies are reproduced by 
assuming that M* is a constant fraction of the halo mass M. The upper bound to the 
low-luminosity LSBs coincides with the virial velocity of haloes in which supernova 
feedback could significantly suppress star formation, V < 100 km s -1 (Dekel & Silk 
1986). We argue that the energy fed to the gas obeys E$n oc M* despite the radia- 
tive losses, and equate it with the binding energy of the gas to obtain M*/M oc V 2 . 
This idealized model provides surprisingly good fits to the scaling relations of low- 
luminosity LSBs and dwarfs, which indicates that supernova feedback had a primary 
role in determining the fundamental line. The apparent lower bound for galaxies at 
V ~ lOkms^ 1 may be due to the cooling barrier at T ~ 10 4 K. Some fraction of the 
dark haloes may show no stars due to complete gas removal either by supernova winds 
from neighboring galaxies or by radiative feedback after cosmological reionization at 
z; on . Radiative feedback may also explain the distinction between dwarf spheroidals 
(dE) and irregulars (dl), where the dEs, typically of V 30kms -1 , form stars before 
z ion and are then cleaned out of gas, while the dls, with V > 30kms _1 , retain gas-rich 
discs with feedback-regulated star formation. 

Key words: galaxies: dwarf — galaxies: formation — galaxies: fundamental param- 
eters — galaxies: local group — supernova remnants — winds, outflows 



1 INTRODUCTION 

The galaxies can be crudely divided into two main classes 
based on their location in the plane of surface brightness 
versus luminosity or stellar mass M* [see for example Fig. 1 
in Dekel & Silk (1986, hereafter DS), Fig. 7a of Kaufhnann 
et al. (2003b, hereafter K03), and other references cited be- 
low]. The bright galaxies, dominated at the bright end by 
ellipticals and early- type spirals, have relatively high surface 
brightnesses which are only weakly correlated with stellar 
mass. The fainter galaxies, spanning the broad range from 
relatively bright late-type spirals all the way down to the 
Local-Group dwarf galaxies, have their conditional average 
surface brightness at a given M* decrease with decreasing 
M,. The transition occurs near M„ ~ 3 x 10 10 M Q (corre- 
sponding to absolute magnitude of about —20.8 and —19.0 



in the r and b bands respectively). Other global galaxy 
properties, such as the mean metallicity, behave in a sim- 
ilar manner as a function of stellar mass, with the tran- 
sition seen at a similar characteristic scale (see references 
below). We refer to these two general classes hereafter as 
HH versus LL galaxies, standing for High-luminosity High- 
surface-brightness (HSB) galaxies versus Low-luminosity 
Low-surface-brightness (LSB) and dwarf galaxies. For the 
purpose of our current idealized theoretical modeling, we 
simply distinguish between these two coarse types based on 
the M» transition scale. 1 This kind of classification can be 

1 LSBs are sometimes defined in the literature by central blue 
surface brightness > 22 mag arcsec - 2 , 1-a from Freeman's (1970) 
mean value for bright HSB spirals (21.65 ± 0.35), and sometimes 
by > 23 mag arcsec -2 . By referring to LL and HH we do not imply 
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traced back, e.g., to Binggeli, Sandage & Tarenghi (1984), 
Wirth & Gallagher (1984), Kormendy (1985) and Hoffman 
et al. (1985). Dekel & Silk (1986) have highlighted this clas- 
sification scheme (their Fig. 1) in the context of their early 
theoretical modeling (their Fig. 6). 

1.1 HH and LL galaxies 

The analysis by K03 of 80,000 galaxies from the Sloan Digi- 
tal Sky Survey (SDSS) highlight the bivariate distribution of 
relatively bright galaxies in the plane of surface brightness 
and stellar mass, above their claimed completeness limits 
of absolute r magnitude —17 and effective surface bright- 
ness 23 mag arcsec -2 . With the spectral information avail- 
able for SDSS galaxies, their stellar masses can be evaluated 
more reliably than before using population synthesis models 
(Kauffmann et al. 2003a). The transition scale shows very 
clearly at M* ~ 3 x 10 10 M© (K03, Figure 7a). The bright 
galaxies in the range 3 x 10 10 < M, < 1O 12 M have their ef- 
fective surface brightnesses scattered about a mean value of 
/i* ~ 1O 9 M0 kpc -2 (referring to the mean surface brightness 
within the half-light radius; the central surface brightness is 
typically larger by a factor of ~ 3), with only a weak sys- 
tematic trend of roughly fi, oc Mf 2 . (When viewed as a 
function of luminosity, the surface brightness is actually de- 
creasing slowly with luminosity in this range, because the 
stellar mass-to-light ratio is increasing, see Blanton et al. 
2003.) On the other hand, the correlation at the top part 
of the LL regime, 10 8 < M, < 3 x 10 10 M Q , is well fit by 
fi* oc Mf 6 (or even slightly steeper). A similar correlation, 
with a slope ~ 0.6 — 0.7, is found in the top LL regime 
between surface brightness and i absolute magnitude based 
on 144,609 SDSS galaxies (Blanton et al. 2003), indicating 
that the translation to stellar mass by K03 makes a neg- 
ligible difference in this regime. A consistent correlation is 
measured from other samples of galaxies as well, e.g., Cross 
et al. (2001) find a slope of ~ 0.42 in the 2dF survey over the 
whole range of galaxies brighter than Mb — —16, de Jong & 
Lacey (2000) measure a slope of 0.5 for Sdm galaxies, Driver 
(1999) finds a slope of 0.67 for a sample of fainter galaxies 
in the Hubble Deep Field, and Ferguson & Binggeli (1994) 
reported a slope of 0.7 for Virgo dwarf galaxies. 

The spread in surface brightness at a given luminosity, 
which is not directly relevant for our theoretical analysis in 
the current paper, is a matter of debate among the observers. 
The uniformly selected SDSS data show a relatively tight 
distribution about the mean relation in the — M* plane, 
both above and below the transition scale (see K03, Figs. 7, 
and 10). We learn for example that low- luminosity galaxies 
with high surface brightness, such as M32, seem to be rare. 
If the photometric completeness limit of the SDSS data in 
r is indeed below 23.0 mag arcsec -2 , then, for a given M* 
near M* ~ 10 9 Mq, this data indicate a significant drop in 
the galaxy count as a function of decreasing surface bright- 
ness. 2 A similar conclusion, of a dearth of high-luminosity 

that there are no LSB galaxies of high luminosity (see below). The 
term "dwarfs" is also used in different ways, referring alternatively 
to small, faint, and low surface brightness galaxies. 
2 Blanton et al. (2003, in preparation) provide further evidence 
for their completeness down to below 23.0 by showing that when 
the galaxies are binned according to Galactic extinction E, the 



galaxies with low surface brightness, is obtained from the 
2dF survey by Cross et al. (2001). On the other hand, there 
are claims in the literature for a significant population of 
low surface brightness galaxies with high luminosities (e.g., 
Disney 1976; Phillipps, Davies & Disney 1988, 1990; Davies 
et al. 1994; Impey et al. 1996; Bothun, Impey & McGaugh 
1997; Sprayberry et al. 1997; papers in Davies, Impey & 
Phillipps 1999; O'Neil & Bothun 2000). 3 The SDSS and 2dF 
results argue that this is either a smaller population, or a 
separate population that occupies a different locus in the 
jU» — M* plane, below the surface-brightness limits of these 
surveys. Independently of this ongoing dispute, we seek in 
the current paper theoretical understanding for the general 
correlation between average surface brightness and stellar 
mass, as indicated by SDSS and the other datasets across 
5 decades in M», from the transition scale to the smallest 
dwarfs. As mentioned in [|8J the scatter in surface brightness 
at a given luminosity may be affected by other physical pro- 
cesses not studied here in detail — the obvious one being 
angular momentum (e.g., Dalcanton, Spergel & Summers 
1997). 

We notice in passing that, given the observed Schechter 
luminosity function and the transition at M* ~ 3 x 1O 1O M0 
(about a factor or 2 below the mass corresponding to 
Schechter's characteristic L*), more than 95% of the galaxies 
are below the transition scale, while most of the light still 
comes from the bright galaxies above this scale. In terms 
of mass, if the mass function of haloes is similar to that 
predicted for the ACDM cosmology by simulations or Press- 
Schechter-like approximations, then the vast majority of the 
virialized mass is in haloes of LL galaxies. This clearly moti- 
vates a major theoretical effort aimed at understanding the 
origin of the low- luminosity LSBs. 

At least two additional independent relations, beyond 
the jU* — M* relation, are apparent in the SDSS data of 
the relatively bright galaxies. The second, based on prelim- 
inary reports (Tremonti et al. 2003, in preparation), is a 
re-confirmation of a scaling relation involving the metallic- 
ity Z, of roughly Z oc M*' s at the high end of the LL regime, 
turning into no significant correlation between Z and M* in 
the HH regime, Z ~const. Similar gradients in gas metal- 
licity have been seen before in other samples of LSBs (e.g., 
McGaugh 1994; de Blok & McGaugh 1997). 

The third, which is typically the tightest correlation 
obeyed by galaxies, is between their luminosity and the char- 
acteristic velocity V measuring the depth of the gravita- 

peak of the apparent surface brightness distribution shifts accord- 
ingly, from about 21.2 at E = to about 23.0 at E = 1.8 (a test 
first applied by Davies et al. 1993 to the ESO-Uppsala catalog). 
3 O'Neil & Bothun argue for a population of luminous LSBs 
based on their finding that the surface brightness distribution 
function is flat down to below a central blue value of 24 - down 
to the survey completeness limit. They infer that many of the 
LSBs are luminous and extended based on the moderate differ- 
ences in the distributions of sizes and velocities between HSBs 
and LSBs in these surveys. With the correlation fi t oc M?' 6-0 ' 7 , 
the expected variation in radius across this decade of [i* is only 
a factor of ~ 2, which may be consistent with the data. The cor- 
responding variation in velocity, V oc M®' 2 , is similarly weak. 
The flat surface-brightness distribution can therefore be consis- 
tent with the correlation detected in SDSS, 2dF and other surveys 
between surface brightness and luminosity. 
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tional potential well, corresponding to roug hly V oc Mf 
(e.g., Bernardi et al. 2003). This is the Tully-Fisher relation 
for the rotation velocity in discs and the Faber-Jackson rela- 
tion (or a projection of the generalized Fundamental-Plane 
relation) for the dispersion velocity in spheroids. We term 
this kind of relation between M, and V a "TF" relation. 
A similar TF relation seems to extend down at least to the 
top part of the LL regime, with no obvious change at the 
transition between HH and LL galaxies (e.g., Zwaan et al. 
1995; Sprayberry et al. 1995; Dale et al. 1999). 

K03 (Fig. 9 and 11) also find that the luminosity "con- 
centration" within the galaxies correlates with stellar mass, 
which they interpret as a measure of bulge-to-disc ratio. The 
bright galaxies are dominated by elliptical galaxies, the LLs 
near M, ~ 10 9 M© are dominated by discs, and in between 
the bulge-to-disc ratio is gradually decreasing with decreas- 
ing M„. Associated with this trend is an increasing gas-to- 
star ratio (e.g., Longmore et al. 1982; McGaugh & de Blok 
1997) and a younger, bluer stellar population in galaxies of 
decreasing M, down to ~ 1O 8 M0 (e.g., de Blok, van der 
Hulst & Bothun 1995; Sprayberry et al. 1995; McGaugh, 
Schombert & Bothun 1995; another interpretation in van 
den Hoek et al. 2000). At a given M, , the galaxies with 
lower bulge-to-disc ratio and younger stellar ages tend to be 
of lower surface brightness (K03, Fig. 14). Preliminary re- 
sults from the SDSS data (Brinchmann et al. 2003, in prepa- 
ration) indicate that while the current star formation rate 
in HHs shows no clear correlation with stellar mass, there is 
a correlation of the sort M 4 oc M, in the top LL regime. 

1.2 Local-Group dwarfs 

The dwarf galaxies of the Local Group (LG, see JSJ, as 
well as those observed in the Virgo cluster and the Local 
Supercluster (e.g., Binggeli & Cameron 1991; Ferguson & 
Binggeli 1994), seem to be in a crude sense an extension of 
the LL population of SDSS, stretching its range to almost 5 
decades in stellar mass: 6x 10 5 < M« < 3x 10 10 M Q . The LG 
dwarfs themselves span almost the whole LL range, up to 
~ 1O 1O M0, and they therefore constitute a very useful sam- 
ple, which we use extensively in this paper. As summarized 
below and displayed in figures in JBJ the correlations between 
the observed global quantities seem to be roughly consis- 
tent throughout this range, including the surface brightness, 
metallicity and velocity as a function of stellar mass. The 
correlations between every pair of parameters for galaxies in 
this regime define a fundamental line in the multi-parameter 
space, say M 4 , /i*, V and Z. For the purpose of our ideal- 
ized modeling of this mean line, we ignore the differences 
between different types within the galaxies that populate 
the line, and any possible population outside this family. 

Within the broad LL family, the brighter dwarfs are 
typically disc-like or irregular in shape and gas rich (we use 
hereafter the broad term "dwarf irregulars", dl); in many 
ways they seem to represent a continuous extension of the 
LLs observed by SDSS. On the other hand, the faint end, 
M» ^ 3 x 1O 7 M0, is dominated by dwarf ellipticals or 
spheroidals (hereafter dE), typically with only little gas and 
current star formation. While the dEs extend the general 
LL trends of decreasing surface brightness and metallicity, 
they clearly represent a reversal in the trend of galaxy type 
(bulge-to-disc ratio, gas-to-stars ratio and star-formation 



rate), and as we show below they seem to show an addi- 
tional difference in the velocity trend (TF). After address- 
ing the origin of the LL/dwarf family as a whole in the main 
body of this paper, we also discuss the possible origin of the 
distinction between dEs and dls. 

The data on the ~ 40 dwarf galaxies of the Local Group, 
based primarily on the compilations by Mateo (1998) and 
van den Bergh (2000), have been analyzed by Woo & Dekel 
(in preparation, hereafter WD). They computed the corre- 
sponding scaling relations in the range 6 x 10 5 < M, < 
1O 1O M0 in comparison with the SDSS results for the bright 
end of the LLs. The stellar mass M, of each galaxy has been 
derived from the observed magnitudes using the mean age 
and metallicity of the stellar population and a simple popu- 
lation synthesis model (kindly provided by G. Kauffmann); 
the results were found to be quite insensitive to the details of 
this derivation. For the central surface brightness WD find a 
tight correlation about the scaling relation ^l, oc 

M 0.55±0.03 j 

extending down to /i, ~ 3 x 1O 6 M0 kpc -2 at the faint end. 
The best-fit slope is determined by WD via a linear regres- 
sion of the log variables, taking into account the errors in 
both (i.e., minimizing the 2-dimensional x 2 as in Numerical 
Recipes, Press et al. 1992, §15.3). The Pearson correlation 
coefficient (NR, eq. 14.5.1) is r = 0.88. The slope is quite 
similar to the slope obtained for the effective surface bright- 
ness of LLs sampled by SDSS, while the amplitude for the 
effective is lower by a factor of ~ 3. 

For the metallicity Z, WD took the stellar [Fe/H] 
(mostly for dE) and/or a constant factor times the Oxygen 
abundance of the gas (mostly for dl), where the constant fac- 
tor has been chosen to minimize the scatter in the Z — M* 
relation. They find a tight correlation, with the best-fit scal- 
ing relation Z oc M°- 40±0 ' 02 and r = 0.92 

For the internal velocity V, WD adopted the observed 
maximum circular velocity for the dls and y/3a p for dEs, 
where cr p is the observed projected central dispersion veloc- 
ity. When the fit is performed across the whole dwarf range, 
the TF scaling relation is V oc M» 24±0 01 with r = 0.89. 
When inspected more carefully, the dls at the bright end 
show a slight steepening which merges smoothly into the 
known TF relation for bright galaxies. At the faint end, 
M t < 3 x 1O 7 M , there is an indication that the veloci- 
ties of the dEs are bound from below by V > lOkrns -1 , 
and can actually be fit by V ~ const, (see [Q. Nevertheless, 
the tight scaling relations over the whole range indicate that 
the LL/dwarf galaxies basically constitute a one-parameter 
family, which calls for a simple physical explanation. 

1.3 Supernova feedback 

The dwarf galaxies are central players in one of the main 
problems facing galaxy-formation theory in the context of 
CDM cosmology — the so called "missing dwarf problem" . 
This refers to the apparent discrepancy between the pre- 
dicted abundances of halo masses in the CDM cosmology, 
especially subhaloes within larger haloes, and the relatively 
few, faint dwarf galaxies observed, e.g., in the Local Group 
(Klypin et al. 1999; Moore et al. 1999; Springel et al. 2001). 
This problem is related to the fact that for galaxies fainter 
than L, the luminosity function is observed to be flatter 
than the mass function predicted for haloes in the ACDM 
cosmology. This implies that the stellar-to-virial mass ra- 



4 A. Dekel & J. Woo 



tio M*/M must be decreasing with decreasing M, namely 
fewer stars were formed per unit total mass in fainter LSB 
galaxies. The systematic variation in M*/M can serve as 
a clue for understanding the LL-LSB phenomenon. Follow- 
ing preliminary ideas by Larson (1974), Dekel & Silk (1986) 
studied the general scenario where the key physical process 
governing this phenomenon is the supernova feedback from a 
first generation of stars, which either drives out a significant 
fraction of the original halo gas or suppresses star formation 
in the retained and added gas. DS showed that the observed 
scaling properties of dwarf galaxies can be qualitatively con- 
sistent with this picture, provided that the potential wells 
are dominated by non-gaseous dark haloes with a structure 
that crudely resembles the predictions of the CDM scenario. 
They studied the amount of energy fed into the interstel- 
lar gas by supernova ejecta subject to radiative losses and 
found that haloes with virial velocities significantly lower 
than a critical value of order ~ 100 km s -1 can lose signifi- 
cant amounts of their gas and/or effectively suppress further 
star formation. 

Recent observations provide cumulative evidence for 
massive outflows from galaxies, consistent with being gen- 
erated by supernovae. For example, HI maps of large disk 
galaxies show empty bubbles associated with outflowing 
clouds, which indicate an energy equivalent to ~ 100 normal 
supernovae (Boomsma et al. 2002). Chandra x-ray measure- 
ments show extended winds of soft x-ray about galaxies such 
as M82, associated with hard x-ray sources in the galaxy, in- 
dicative of young massive stars and supernovae (Roy et al. 
2000; Martin, Kobulnicky & Heckman 2002). Outflows are 
detected directly in local starburst galaxies, and are seen to 
be driven by SNII activity (e.g., Legrand et al. 1997; Mar- 
tin 1999b; Heckman et al. 2001). In some of these galax- 
ies, outflows of a few hundred kms -1 are inferred based on 
the blueshifted metal absorption lines from the approach- 
ing foreground compared to the redshifted Ly-a photons 
backscattered from the receding background. Spectroscopy 
of the brighter lensed galaxies at high redshifts reveal using 
a similar effect typical outflows of 200 — 800 km s" 1 (Franx 
et al. 1997; Frye & Broadhurst 1998; Pettini et al. 2001; 
Frye, Broadhurst & Benitez 2002). New evidence for strong 
outflows at z ~ 3 is provided by Adelberger et al. (2003), 
who interpret their measurements of absorption systems in 
QSO spectra as bubbles of radius ~ l/i _1 Mpc (comoving) 
around Lyman-break galaxies, almost empty of neutral hy- 
drogen. Combined with the measured outflow velocities of a 
few hundred kms -1 , this is an indication for energetic winds 
that persist for a few hundred million years and drive away 
nearby intergalactic gas. A strong observed correlation of 
galaxies with intergalactic metals supports the idea that the 
IGM has been enriched by the outflows from Lyman-break 
galaxies. In several cases, the outflow rate is observed to be 
proportional to the star formation rate (Martin 1999a), con- 
sistent with a stellar feedback origin for the outflows. This 
body of evidence indicates that supernova-driven winds ac- 
tually exist, which helps motivating our theoretical modeling 
of the relevant features of galaxy formation. 

In this paper we improve the DS scenario for the for- 
mation of dwarf galaxies in view of the developments in 
galaxy formation theory and the refined observed scaling re- 
lations across the LL/dwarf family. Using a simple energetics 
criterion and standard assumptions regarding the origin of 



galaxy sizes, we now show that the observed scaling relations 
naturally emerge from the simplest possible supernova feed- 
back scenario, even before one tries to model and simulate 
in detail the complex physics of the feedback mechanism, 
and before one worries about the different types of dwarf 
and LSB galaxies and the scatter in their properties. We 
then address the possible role of radiative feedback in dis- 
tinguishing between dE and dl galaxies, and in preventing 
star formation altogether in some halos. 

In J2]we address the role of standard assumptions in de- 
termining the scaling relations for galaxies in general. The 
assumptions include spherical collapse to virial equilibrium 
in CDM haloes and angular-momentum conservation. In 
we apply the analysis to bright galaxies where we assume 
that M,/M ~ const.. In 0we summarize the DS deriva- 
tion of the velocity characterizing the supernova-feedback 
scale. In JHJwe use simple theoretical considerations regard- 
ing supernova feedback to derive the scaling relations of the 
LL/dwarf family. In ^we compare the model predictions 
to the observed relations shown by the Local Group dwarfs, 
and comment on the comparison with the SDSS data. In J7| 
we discuss the possible role of radiative feedback in distin- 
guishing between dEs and dls. In JUwe discuss our results 
and related issues. 



2 SCALING RELATIONS: GENERAL 

We show that the basic observed scaling relations for galax- 
ies in the two regimes can be reproduced to a surprising 
accuracy based on the simplest possible physical assump- 
tions. These include the virial theorem for spherical cold- 
dark-matter haloes, and the notion that a fraction r\ of the 
original gas makes stars in a disc such that the size of the 
stellar system is determined by angular momentum. For HH 
galaxies we recover the scaling relations by taking n to be 
independent of halo mass, assuming that feedback is not too 
effective there. For LL galaxies, where feedback is a key fac- 
tor, we use in SJ3 below a simple energy constraint for the 
efficiency of supernova feedback to predict how rj should vary 
with the halo virial velocity. Together with the constraints 
from the virial theorem and angular momentum, this leads 
to the characteristic scaling relations in the LL regime. In 
this section we derive the scaling relations in general terms 
without specifying the behavior of rj. 

Assume a dark-matter halo of mass M reaching virial 
equilibrium at a time corresponding to cosmological expan- 
sion factor a = (l + z) -1 . The virial radius R is defined in the 
spirit of the spherical collapse model by a given density con- 
trast A relative to the mean universal density at that time, 
namely by M/R 3 oc A a -3 . At early times, when f2 m ~ 1, 
the relevant density contrast is A ~ 180, while for the stan- 
dard ACDM cosmology (with £7a = 0.7 and fi m = 0.3 to- 
day) it rises to A ~ 340 today. In the following, we ignore 
the weak redshift dependence of the A factor. 4 The virial 

4 The maximum change is obtained at low redshifts. For example, 
in the range 2 = — 2 the change is roughly A oc a 1 / 2 , which 
implies that a in the following expressions should be replaced by 
A~ 1//3 a oc a 5 / 6 . This is a weak effect, which becomes even weaker 
at higher redshifts. 
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velocity is defined by the virial theorem, V 2 oc M/R, such 
that the three virial quantities at a define a one-parameter 
family: 



M oc a 3/2 V 3 oc a~ 3 R 3 . 



(1) 



In the simplest analysis we ignore the possible system- 
atic increase of a as a function of halo mass (but see below) . 
This dependence is relatively weak already as predicted by 
cosmological spherical collapse in the ACDM cosmology, it 
gets weaker for smaller haloes as the rms density fluctua- 
tions approach a constant on small scales, and it is weaken 
further by effects like the merging of early-forming small 
haloes into bigger ones (see Wechsler et al. 2002). The virial 
relations for typical haloes thus take in this approximation 
the simple form 



M oc V 3 oc R 3 . 



(2) 



We may keep tracing the a dependence in the general ex- 
pressions below in order to allow small corrections due to its 
possible weak dependence on M, when desired. 

Considering next the baryonic component, we assume 
that the halo is initially filled with gas of mass M g ~ fbM, 
where /b (— 0.13) is the universal baryonic fraction. For 
large galaxies, the gas is assumed to be shock-heated to the 
halo virial temperature, but as long as M < 10 12 — 1O 13 M0 
the gas in the halo can cool in a dynamical time (shorter 
than the Hubble time) and contract to form stars (Rees & 
Ostriker 1977; Silk 1977; White & Rees 1978; Blumenthal et 
al. 1984 in the context of dark haloes). We denote the ratio 
of luminous stellar mass M* to initial gas mass M s by 77, 



M» = rjM g oc 7]M, 



(3) 



without yet specifying how 77 may depend on M. Substi- 
tuting eq. in the virial relations, eq. Q, we straightfor- 
wardly obtain a general TF relation between V and M*: 



-1/2 -1/3,^1/3 
V oc a 1 77 ' M, ' . 



(4) 



As long as the halo rotation curves are roughly flat at large 
radii, we ignore the difference between the virial velocity V 
and the observed velocity Vmax- 

If the baryons within the halo virial radius R cool and 
contract to a centrifugally supported disc of radius R* while 
preserving their specific angular momentum j, then, follow- 
ing Fall & Efstathiou (1980) and Mo, Mao & White (1998), 
we write i?* ~ XR, where A = j/ (RV) is the initial baryonic 
spin parameter (according to the revised, practical defini- 
tion of Bullock et al. 2001b). Then, from the virial relations 
above, 



1/3 



R* oc XaM 

With M* oc r\M, this implies for the surface brightness 
fi* oc M t R~ 2 oc A" 2 a" 2 T/ 2/3 M, 1/3 . 



(•») 



(6) 



The characteristic radii and surface brightnesses derived 
for discs can be argued to be roughly valid also for the 
spheroidal stellar components, elliptical galaxies or bulges 
of spirals. This is based on energy conservation and the 
virial theorem under the assumption that the spheroids were 
formed by mergers of discs. 

The distribution of halo spin parameter is known from 
cosmological simulations to be insensitive to halo mass 
(Barnes & Efstathiou 1987; see Bullock et al. 2001b for 



the case of ACDM cosmology), so for a crude approxima- 
tion it could be dropped from the mean scaling relations 
for the radius and surface brightness (but see a comment 
in S0. However, we note in passing that the A dependence 
in the expressions for the radius and surface brightness, as 
opposed to its absence from the TF relation, may have an 
important implication on the scatter about these relations. 
Cosmological simulations of ACDM reveal that the distri- 
bution of halo spin parameter is log-normal with a mean of 
A = 0.035 ± 0.005 and a standard deviation in the (deci- 
mal) log of cr A = 0.184 ± 0.011 (Bullock et al. 2001b). This 
is already comparable to the standard deviation of the dis- 
tribution of stellar radii for the LL galaxies in SDSS (K03, 
Fig. 8), even before considering the additional scatter in a 
and in 77. The absence of A dependence in the TF relation, 
and its weaker dependence on a and 77, may explain why the 
TF relation is much tighter (see Courteau et al. 2003). 

The amount of metals produced in a galaxy is assumed 
to be proportional to M* with a constant yield y. When 
r\ ~ 1, the metallicity is simply Z ~ y. When 77 <C 1, in 
the instantaneous-recycling approximation (e.g., Searle & 
Sargent 1972; Audouze & Tinsley 1976), one expects Z ~ 
y-q. So we approximate in general 



Z oc M«/M g 



(7) 



As mentioned above, small corrections to the above re- 
lations may result from a correlation between the time of 
formation of a halo and its mass. An upper bound to this 
effect may be obtained from the straightforward prediction 
based on spherical collapse, ignoring the fact that many 
early-forming small haloes eventually merge to bigger ones 
and thus weaken the a(M) relation. For a power spectrum 
of linear density fluctuations that resembles the power law 
Pfe oc k n at the vicinity of the scales relevant for galactic 
haloes, the typical mean density fluctuation within a proto- 
halo is 5 oc M~ {n+3)/(i D(t), where D(t) is the linear growth 
rate, D(t) oc a for the Einstein-deSitter cosmology relevant 
at high redshifts. The formation time in the spherical col- 
lapse model can be approximated by S ~ 1.7 for the linearly 
extrapolated mean density fluctuation, so one obtains 



a oc M 



(n+3)/6 



(8) 



The virial relations for typical haloes, eq. ifT). thus become 
The TF relation, eq. @, is now 



V oc (y~ M„) 



1 m ^i-™)/ 12 



(10) 



The stellar radius is now given by R, oc \M f5+n ^ 6 such 
that the surface brightness, eq. @, is replaced by 



H, oc A-y 5+n)/3 M- (2+n)/3 . 



(11) 



The expressions so far should be valid in general, both 
for HH and LL galaxies. The differences between the two 
classes enter mainly via the behavior of 77, with an additional 
weak effect due to the difference in the effective n in the 
maximum M dependence of a. 
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3 HH GALAXIES 



£ S N(i) evM,t lad f(t), 



(14) 



For HH galaxies we take r\ to be roughly independent of 
halo mass. This is based on the assumption that feedback 
effects do not significantly heat or remove most of the gas 
from these galaxies such that most of the gas, or a constant 
fraction of it, eventually forms stars (see JIJ. With r/ and 
a independent of mass in eq. (|1J, eq. ijfjj and eq. 0, the 
scaling relations for HH galaxies become 



VocMt 



1/3 



jit* oc M* 



1/3 



const. 



(12) 



These are already in qualitative agreement with the observed 
relations for HH galaxies. 

When considering the limit of maximum M dependence 
of a, we recall that big galactic haloes in a ACDM cosmology 
correspond to the part of the power spectrum where n ,$ — 
2. For example, with n = —2 at the bright end, one has 
a oc Af 1 / 6 . Then the virial relations become M oc V 4 oc 
R 2 . With r\ assumed independent of V for HH galaxies, and 
with the maximum M dependence of a computed above, the 
predicted scaling relations become 



Vmax OC M, 



1/4 



const. 



const. 



(13) 



In order to compare with observations in terms of luminosity 
rather than stellar mass, one can assume that for HH galax- 
ies the stellar mass-to- light ratio varies like M*/L oc L 0,3 
(e.g., Courteau et al. 2003). [This is based, for example, 
on the reading of Fig. 7 of Courteau et al. that (V — I) ~ 
-0.09M/, combined with the result from Table 1 of Bell & 
de Jong (2001) that log(M*/Li) ~ l.35(V - I). A similar 
result is obtained for elliptical galaxies (Bender, Burstein & 
Faber 1992).] With Vmax — V we thus roughly recover the 
observed TF relation: L oc V^^. The surface brightness as 
measured in terms of luminosity is predicted to be slowly 
decreasing with luminosity: / oc L/R 2 oc L -0,3 , in qual- 
itative agreement with observations. We note that in the 
TF relation the correction due to the correlation of a and 
M roughly balances the correction due to the correlation of 
M*/L and M. However, the corresponding corrections to the 
relation of surface brightness and luminosity add up. In any 
case, the corrections to the simple predictions of eq. 11211 are 
small, and we expect the predicted scaling relations for HH 
galaxies to lie somewhere between the relations in eq. 11 -t 
and eq. 11131 . This range is in general agreement with the 
observed scaling relations for HH galaxies. 

Recall that, beyond the standard assumptions of virial 
equilibrium and spherical collapse, the key assumption for 
HHs was that most of the original gas, or a constant fraction 
of it, turns into stars, namely, r\ ~ const. 



4 SUPERNOVA FEEDBACK: THE CRITICAL 
SCALE 

Dekel & Silk (1986) evaluated the maximum total energy 
fed into the interstellar gas by a collection of supernova ex- 
plosions due to a period of star formation at a constant 
rate M», taking into account the radiative loses (based on 
the standard evolution of a supernova remnant in a uniform 
interstellar medium, e.g. Spitzer 1978; Ostriker & McKee 
1988). They found that at time i this energy can be approx- 
imated by 



where e is the initial energy released by a typical supernova 
(e ~ 10 51 erg) and v is the number of supernovae per unit 
mass of forming stars (which for a typical IMF is v ~ 1 per 
50Mq of stars). The characteristic time i ra d marks the end of 
the "adiabatic" phase and the onset of the "radiative" phase 
of a typical supernova remnant, by which it has radiated 
away a significant fraction of its energy. The dimensionless 
factor f(t) turns out to be of order unity when i ~ i ra d; it 
grows roughly oc i for t < i ra d and oc t 0A for t > i ra d- 

We assume here that the stellar population of mass M* 
has formed over some constant multiple r of the free-fall 
time tg, namely 



M, 



(15) 



Substituting in eq. (I14H we obtain that the total energy fed 
into the gas is 



i? SN ocAf»^. 

iff 



(16) 



DS noticed that in the temperature range 6 x 10 4 < T < 
6 x 10 K the cooling rate scales approximately like A oc T _1 , 
which implies that the ratio i ra d/te is roughly a constant, 
of order 10 -2 , independent of the gas density or the halo 
parameters. This leads to Esn oc M*, which we show below 
is a key for deriving the scaling relations of LL galaxies (|HJ. 
Note that DS originally assumed M* oc Mg/tg rather than 
the oc M t /tg of eq. (11511 . The two assumptions are roughly 
equivalent in the case of bright galaxies and when trying to 
estimate the transition scale between HH and LL galaxies 
where M* ~ M s . 

DS also showed that if star formation is rapid, r ~ 1, 
then the filling factor of the expanding supernova shells 
within the halo is of order unity when the typical shell is 
at the end of its adiabatic phase, at t ra d. This coincidence 
indicates that the supernova energy (minus the radiative 
losses) can be fed quite evenly and efficiently into most of 
the gas via the expanding shells that reach a significant mu- 
tual overlap roughly at the time after which they become 
ineffective. It also justifies the adoption of / ~ 1 in eq. 1141 . 

A necessary condition for heating or unbinding most 
of the initial gas of mass M s is obtained by requiring that 
the energy fed by supernovae is comparable to the binding 
energy of the gas in the halo potential well, 



E S n = (l/2)M g V 2 



(17) 



Here V is the virial velocity of the halo, which we assume for 
simplicity to be isothermal and to dominate the potential. 
DS then pushed this approximate relation to the limit where 
a large fraction of the gas turns into stars and obtained the 
critical velocity 



Vsn 



„ „ V t ra d 
^J e 7 

t iff 



1/2 



100 kms" 



(18) 



This critical velocity is evaluated using the typical values of 
e and v with / ~ r ~ 1 and is independent of the gas density 
because of the robustness of i ra d/iff- The interpretation of 
this critical velocity is that gas removal becomes possible in 
haloes with virial velocities smaller than Vsn- We note that 
the corresponding virial mass is 
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M SN ~ 2.2 x lO n M ( 



Vsn 



100 km s" 1 



,2/3 



(19) 



With M* ~ M g ~ fhM and the universal baryonic frac- 
tion /b — 0.13 the corresponding characteristic stellar mass 
today is 



M, SN ~3x 10 1U M Q , 



(20) 



in excellent agreement with the transition scale seen in the 
SDSS data. 

The haloes of bright galaxies, which have retained most 
of their gas, are thus limited to the regime of deep poten- 
tial wells, V > Vsn. We associate the galaxies that form in 
haloes below the critical supernova scale with LL or dwarf 
galaxies. The importance of feedback effects in the history of 
these galaxies implies that their scaling relations can be very 
different from those shown by galaxies that live in haloes of 
virial velocities larger than Vsn- 



5 LL GALAXIES 

We now use the feedback energetics constraints, eq. 1161 and 
eq. 11711 . to determine the behavior of 77 in the LL regime, 
where we expect feedback effects to allow only a fraction of 
the gas to turn into stars, 



77 ; 



< 1. 



(21) 



The supernovae resulting from the first burst of stars either 
blow out the rest of the gas, or at least provide enough 
feedback energy to regulate the subsequent star formation 
rate and keep it low. We assume that 77 ~ 1 at V = Vsn ~ 
100 km s -1 and that 77 becomes gradually smaller for haloes 
of smaller velocities. The following simple analysis is actually 
carried out in the limit of strong feedback, 

Our key starting point is eq. 1161 with trad/iff = const., 
namely the energy fed into the interstellar gas by supernovae 
is proportional to the final stellar mass, 



E SN oc M, . 



(22) 



Without the radiative losses of the supernova energy, this 
would have been anybody's first intuitive guess for a relation 
between these quantities. We argue here, in the spirit of 
the DS analysis, that the actual energy fed into the gas 
after significant radiative losses is still a constant fraction of 
the original supernova energy. This makes eq. 1221 a valid 
approximation in the realistic case. 

In order to allow significant heating or total blowout of 
the initial gas, the total input by supernovae should be at 
least comparable to the binding energy of the gas, eq. 1171 . 
With cq. 1221 . the energy condition becomes M» oc M g V 2 , 
namely 



Tj oc V 



(23) 



The scaling relations for LLs all follow from this basic re- 
lation, which measures the strength of the feedback effects 
along the halo sequence characterized by the parameter V 
in the range V < Vsn. 

Eq. 1231 . combined with the virial relations for the halo, 
eq. Q, and then M* = T]M, yield 



Recall that in the instantaneous-recycling approximation, 
for 77 <C 1, the metallicity is simply 



Z oc 77, 



(25) 



77 oc a M 



2/3 



-3/5 pj-2/5 



(24) 



so the mean scaling relation involving metallicity is given by 
eq. l ISH . 

Substituting 77 in eq. Q we obtain for the TF relation 
in the LL regime 

V oc a- 3/10 Mt /s . (26) 

Then substituting 77 in eq. (jy and eq. © we obtain for 
the radius R* oc Aa 6//5 A/* / ' 5 and for the surface brightness 

/i* oc A~ 2 a~ 12/5 M 3/5 - (27) 

In summary, when ignoring possible weak systematic 
dependences of a and A on M, the scaling relations for 
LL/dwarf galaxies are predicted to be 

VocM„ 1/5 , ZocM 2/5 , fi«(xM* /5 . (28) 

In order to evaluate the maximum correction due to the 

possible dependence of a on M, we use eq. (JSJl in eq. 1241 
and obtain 

77 oc M (1 " n)/6 oc Mi l - n)K7 - n) . (29) 
Then the TF relation, eq. 1261 . becomes 

V oc a- 3/10 M^ 5 oc M (i-«)/C"-2n) _ (3Q) 
and the surface brightness, eq. (1271 . becomes 
^ocA- 2 M t - 3(1+n)/(7 " n) . (31) 

Very small galaxies in the ACDM cosmology correspond to 
the part of the power spectrum where n is not much larger 
than the lower limit of n = —3, implying a similar formation 
time for dwarf galaxies of all masses and therefore a constant 
a in the above relations, thus leading to eq. 1281 . For LL 
galaxies not much below Vsn we may try for example a typ- 
ical n = -2.5, for which a oc Mi" +3)/(7 " n) oc M, 1/19 . This 
implies negligible effects on the TF relation and the metal- 
licity relation, but the weak correction to the surface bright- 
ness relation may be marginally detectable, oc A~ 2 M* / ' 19 
compared to (i, oc A _2 M« //5 . 

In this case, however, we may also wish to incorporate 
the possible mass dependence of A. To a first approxima- 
tion, as said above based on cosmological simulations, the 
distribution of halo spin parameter is independent of the 
halo virial properties and its formation time. As long as the 
baryons initially trace the spatial distribution and kinemat- 
ics of the halo, their A distribution can be assumed inde- 
pendent of M and a. However, while the baryons in bright 
disc galaxies seem to have spin parameters similar to those 
of their host haloes, LSB disc galaxies may tend to be as- 
sociated with a higher spin parameter. For example, van 
den Bosch, Burkert & Swaters (2001, hereafter BBS) stud- 
ied the spin in a sample of 14 LSB discs with an estimated 
average of V ~ 60kms . They found an average spin pa- 
rameter about 50% larger than that of the dark haloes (see 
Mailer & Dekel 2002, Fig. 8). At the same time, BBS es- 
timated in these galaxies an average baryonic fraction of 
only /d — 0.035, which translates in our terminology to 
?7(V = 60) = / d //b ^ 0.27. Following Mailer & Dekel 
(2002), we model these systematic trends based on preferen- 
tial blowout of low-spin material in dwarf galaxies. In order 
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to obtain the 50% change in spin parameter over the same 
r] range (between 1 and 0.27) the effect of blowout should 
roughly scale like A oc 77" 03 oc a ' 18 AC ' 12 . Plugged into 
eq. 1271 . using n — —2.5, we now obtain fj,, oc Mj' 7 . This 
kind of correction to the surface-brightness relation should 
be valid for relatively large LL galaxies, where n is not too 
close to —3 and where the BBS analysis indicates a system- 
atic spin dependence. For smaller dwarfs the actual relation 
may be better approximated by eq. 11271 with constant a and 
A, namely (i, oc A'/»' 6 . 

The scatter about the mean scaling relations is expected 
to partly reflect the random scatter about the mean a and 
A. Based on eq. 1271 . the scatter about the mean relation 
H*(M*) is expected to be significant, dominated by the scat- 
ter in A, while the scatter in V and Z is expected to be 
smaller. The residuals in these different relations are ex- 
pected to be correlated. For a given M*, galaxies that lie 
at the bottom of the (j,* distribution are expected to be of 
relatively high a (late formation time) and high A. In turn, 
based on eq. 1261 and eq. 1251 . these galaxies compared to 
the average for that M* are expected to be of low V (though 
high M, given the high a) and low Z. 



6 MODEL VS. LOCAL-GROUP DWARFS 

The success of the simple feedback model for LL galaxies as 
described in J)]can be evaluated by comparing the predicted 
scaling relations, eq. 1281 . to the observed scaling relations 
for LLs in SDSS and in the Local Group. Given the idealized 
nature of the straightforward model, one might only hope 
for a crude qualitative fit. 

The match of the predicted characteristic scale for su- 
pernova feedback, V ~ 100 km s~ 1 , with the observed transi- 
tion at M* ~ 3 x 1O 1O M0 is already remarkable; it indicates 
that this transition may indeed be associated with the onset 
of supernova feedback effects. 

Fig. shows the central surface brightness /Lt* versus 
stellar mass M* for the Local Group dwarfs (from WD). 
The galaxies are either of the two major types, dl and dE, or 
transition cases marked Tr. The data are fit very well by the 
predicted scaling relation /Lt* oc M, ' 6 throughout the whole 
LL range, spanning 5 decades in M*. We do not attempt 
to normalize the predicted relation, and therefore the model 
line in the figure is normalized artificially to provide the best 
fit for the predicted slope of 0.6. The correlation is relatively 
tight, with a Pearson's correlation coefficient for the logs of 
r = 0.88. The model slope is also a good fit to the SDSS 
data in the LL range (Kauffmann et al. 2003b); even the 
predicted slight steepening to oc M*' 7 or so can be seen 
at the bright end of the LL range. The SDSS data refers 
to the surface brightness within the half-light radius, which, 
for an exponential profile, is a factor of ~ 3 smaller than the 
central value. With this relative normalization, the bright 
end of the Local-Group dwarfs lies along the upper 68% 
contour of the SDSS distribution (Fig. 7a of Kauffmann et 
al. 2003b). 

Fig. [5] shows the metallicity Z versus stellar mass M t 
for the Local Group dwarfs (WD). The predicted metallicity 
relation, Z oc Mj' 4 , is a very good fit to the Local Group 
dwarfs. The correlation is tight, with a correlation coefficient 
for the logs of r = 0.92. The preliminary SDSS data indicate 



Figure 1. Central surface brightness versus stellar mass for the 
Local Group dwarfs (from WD). Shown are the regression line 
(1* oc M2' 55 (solid), the correlation coefficient r, and the toy- 
model theoretical prediction fj,, oc M*' 6 normalized for best fit 
(dashed). 



Figure 2. Metallicity versus stellar mass for the Local Group 
dwarfs (WD). Shown are the regression line Z oc M°- 40 and the 
toy-model theoretical prediction Z oc M® A (normalized for best 
fit). 



a similar and perhaps slightly steeper relation at the bright 
end, Z oc M°- 5 . 

Fig- El shows the velocity V versus stellar mass M* for 
the Local Group dwarfs (WD). The predicted relation, V oc 
M*' 2 , is an acceptable eye-ball fit to the data, despite the 
fact that the formal regression slope is somewhat steeper, 
0.24 ± 0.01. The correlation is tight, with r = 0.89. 

We see that the idealized theory for supernova feedback 
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Figure 3. Velocity versus stellar mass for the Local Group dwarfs 
(WD). Shown are the regression line V oc Mj? 24 over the whole 
range and the toy-model theoretical prediction V oc M^' 2 . Note 
the lower bound at V ~ lOkms -1 for the dEs below 3 X 10 7 Mq. 



Figure 4. Surface brightness versus "spin parameter" in three 
bins of constant M, values. The "spin parameter" is actually the 
ratio of stellar radius to halo radius, A cx Rt,/R. The decreasing 
trend is to be compared to the predicted /j* oc A -2 , eq. l6l . 



provides a surprisingly good fit to the characteristic scale 
and to the three independent scaling relations valid across 
the whole LL range. This indicates that the supernova feed- 
back effects, via the parameter r/, indeed have a primary role 
in determining the gross features of the galaxy properties in 
the LL regime. 

The data from SDSS also allow a quantitative evalua- 
tion of the distribution of galaxies about the mean relation 
in the /i»-M, plane (Kauffmann et al. 2003b). In the LL 
regime, the spread in fi, at a given M* is roughly consistent 
with the spread in spin parameter A for haloes of a given 
mass as measured in N-body simulations of ACDM. This 
indicates that A can indeed serve as the main secondary pa- 
rameter for the LL family. 

The idealized theory predicts a dependence of /i* oc 
A~ 2 for a given M* [eq. In Fig. 0] we test the self- 

consistency of this prediction for the Local Group dwarfs. 
We display /i* versus an estimated "spin parameter" given 
by A oc R*/R, the ratio of stellar radius to halo radius. 
The stellar radius is determined from M* and the central 
surface brightness jit*, R* oc (A/*//i*) 1//2 . The halo radius is 
the virial radius corresponding to virial velocity V, where V 
is max{V c irc, x/3o"p} as described in WD. When dividing the 
galaxies into three relatively narrow bins of M» values, we 
see that there is indeed a systematic trend within each bin, 
though slightly flatter than the expected (i* oc A -2 . 

Although the surface brightnesses of dEs and dls follow 
in general a similar scaling relation, the dls do tend to lie 
somewhat below the best fit line. This is consistent with the 
finding in SDSS that at a fixed M, the galaxies with lower 
bulge-to-disc ratio and younger stellar populations tend to 
have a lower surface brightness. These trends are qualita- 
tively consistent with the a and A dependences predicted in 
eq. Cfl . 

While the predicted V oc M°' 2 is a good fit across the 



whole dwarf range, a more detailed investigation of Fig. [£] 
reveals very interesting secondary features. First, there is 
an apparent lower bound for galaxies at V ~ lOkms -1 . 
Second, there is an apparent transition at M, ~ 3 x 1O 7 M0. 
The fainter galaxies can actually be well fit by V ~ const.. 
The dwarfs brighter than 3 x 1O 7 M0 are then fit by a line 
which could be as steep as V oc Mj' 4 . Since these velocities 
are measured in the inner regions of the haloes, they can be 
regarded as lower bounds to the actual dispersion velocities 
of the haloes. If the velocities of the dE haloes are actually 
larger than these lower estimates (as argued by Stoehr et 
al. 2002), then the difference between the TF relation in the 
two regimes, below and above 3 x 10 7 M Q , could become even 
more significant. We note that the faint part is dominated 
by dEs while the brighter part is mostly dls. These are clues 
for the origin of the distinction between these two types of 
dwarf galaxies, which we address in the following section. 



7 RADIATIVE FEEDBACK: dE VS. dl 

After demonstrating the encouraging success of supernova 
feedback in explaining the basic systematic trends in the LL 
family as a whole, we now attempt to consider the possible 
role of another feedback mechanism, and in particular how it 
may differentiate between dEs and dls within the LL family. 

7.1 Radiative feedback 

Cosmological reionization of Hydrogen is complete by Zi on ~ 
6—7 (see a review by Barkana & Loeb 2001, hereafter BL; 
also Loeb & Barkana 2001). The flux of UV radiation that is 
generated by the first stars or AGNs heats and photoionizes 
the gas in the IGM and in virialized haloes (except perhaps 
for the inner regions which can become shielded). As long 
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as the ionizing flux persists, the gas is kept at a fixed tem- 
perature of Tion ~ (1 — 2) x 10 4 K. This can be regarded 
as another feedback mechanism; it can suppress star for- 
mation and clean haloes from gas in two ways. First, by 
photo-evaporation of gas already in haloes. Barkana & Loeb 
(1999) estimated that haloes of V < lOkms" 1 would lose 
most of their gas. Their analysis provides an estimate of 
the gas loss during the first dynamical time after Zj on . How- 
ever, if the gas is kept ionized until z ~ 1 — 2, a dynami- 
cal calculation of evaporation by a continuous wind reveals 
that photo-evaporation would remove the gas from some- 
what larger haloes, up to Vevap — 20 kms" 1 (Shaviv & Dekel 
2003). Second, based on simulations and computations of 
Jeans mass, the pressure of the hot IGM shuts off gas in- 
fall into even more massive haloes, those with velocities up 
to Vjeans — 30kms _1 (see BL §6.5 and references therein). 
Gas could resume falling into small haloes after z ~ 1 — 2 
when the UV background flux declined sufficiently (Babul 
& Rees 1992), but only haloes of V > 20 - 25kms _1 can 
form molecular hydrogen by z ~ 1 and then cool further to 
make stars (Kepner et al. 1997). 

In the presence of a halo potential well characterized 
by a velocity V, the fraction of gas of temperature Tj on that 
is bound to the halo can be estimated by the Boltzmann 
distribution, 

C 1 (-V 2 /kT:„„) f 0n \ 

/bound ocl-e^ ' lon >. (32) 

The velocity corresponding to Ti on is on the order of Vevap ~ 
20kms~ 1 mentioned above. Note that in the limit V < 
Vevap, eq. H32H predicts /bound oc V 2 (as pointed out by 
J. P. Ostriker in a private communication). This reminds us 
of the energy relation for supernovae, 77 oc V 2 , which led 
to the global scaling relations of LLs in J^J Could radia- 
tive feedback (rather than supernova feedback) be the ac- 
tual mechanism responsible for the global scaling relations 
of LLs? First, it is unlikely that a mechanism whose char- 
acteristic scale is ~ 20 — 30kms _1 can be dominant in de- 
termining the observed critical scale of ~ lOOkrns -1 and 
the properties of bright LLs not much below this scale. Sec- 
ond, /bound in the radiative case refers to the sum of bound 
mass in stars and in gas while r\ in the supernova case refers 
to the stellar mass only. While M,/M oc V 2 is consistent 
with the observed scaling relations for LLs, a similar relation 
for the gas-to-mass ratio does not seem to be in agreement 
with the observed trend, especially not in the large, gas-rich 
LLs. There are indications that as one moves from bright to 
fainter galaxies the ratio of gas to stellar mass increases until 
it reaches a maximum at some intermediate scale typical to 
dls before it starts decreasing towards the dE regime (e.g., 
McGaugh & de Blok 1997) . We interpret this as another ev- 
idence against radiative feedback being the dominant mech- 
anism in determining the global properties of galaxies in 
the upper LL regime. The radiative feedback should have 
an important effect though in the small, gas-poor dEs, and 
possibly a complementary effect to the supernova feedback 
in the larger, gas-rich dls and LLs. 

In haloes of V < Vevap, stars can form only before the 
reionization epoch (and possibly much later, at z < 1 — 2). 

If Vevap 

< V < Vjeans, gas that cooled and collapsed before 
the reionization epoch can turn into stars in a slow rate 
also at later times, but new gas cannot be accreted. These 
effects could lead to the gas-poor dEs. On the other hand, 



Figure 5. Velocity versus stellar mass for the Local Group dwarfs 
(WD), same data as in Fig. [3] Shown here is the best fit horizontal 
segment below M, = 3 X 10 7 A/q and the best fit line at larger 
stellar masses: V oc Mj?' 37 . 

in haloes of V > Vjeans there is no much radiative gas loss. 
Galaxies that form in such haloes can retain some gas that 
has not been blown away by supernova winds, or has come 
back after such blowout, and thus give rise to gas-rich dls. 
We thus propose that the main role of radiative feedback is 
to clean up the dEs from their gas and to help regulating 
star formation in dls. 



7.2 Dwarf elliptical galaxies 

Fig-Elshows the same data as in Fig. 01 but with separate fits 
below and above A/* = 3 x 10 7 M©, in the ranges dominated 
by the dwarf spheroidals of the Local Group and by dls 
respectively. In the low-M* range we fit a horizontal line, 

V = 15kms _1 , and then determine the best-fit slope in 
the high-M* range, V oc M° 37 . This fit with a broken line 
naturally provides a better fit than with a single line across 
the whole LL range. It does provide a crude hint for different 
TF relations in the two regimes. 

We thus propose that haloes in the range 10 < V < 
30kms _1 tend to form gas-poor dEs. Efficient cooling by 
Hydrogen recombination at T ^ 10 4 K leads to an early burst 
of stars. The associated supernovae blow out much of the 
gas. The rest of the gas photo-evaporates (if V < Vevap) and 
no new gas can fall in (if V < Vjeans), leaving behind a gas- 
poor system with no significant recent star formation. The 
gas in haloes of V < lOkms -1 cannot cool to form stars at 
any early epoch. This confines all the dEs to a narrow range 
of halo velocities, V ~ 20 ± lOkrns -1 . The spread of M* 
within the dE family is thus predicted not to correlate with 

V (as it does for the dls) but rater to represent variations in 
other quantities such as the time available for star formation 
between the halo collapse and Zi on (as suggested by Miralda- 
Escude & Rees 1998). The spread in M* for a given halo V 
is large for the dEs and much smaller for the dls and the 



Feedback in dwarf/LSB galaxies 11 



spin (see Bullock el al. 2001b). The observed relation of 
roughly fi* oc Af°' 6 tells us that the required trend should 
roue hly be R* oc Xa 1/2 oc M° 2 . If a is indeed anticorrelated 
with M„ for dEs, we expect for A a stronger dependence than 
A oc A/°' 2 . Such a spin gradient could also explain why the 
dwarfs at the faint end are low-spin spheroidals while the 
brighter dwarfs tend to be centrifugally supported discs. As 
a consistency check, Fig.|S|shows the stellar radius 7?, versus 
stellar mass M* for the dEs in the range Af.^x 1O 7 M . 
There is indeed an apparent trend, best fit by R„ oc Af*' 3 
and reasonably consistent with the required i?* oc Xa 1 ^ 2 oc 
M°- 2 . Again, given the limited width of the dE range, we can 
probably tolerate there a certain deviation from the global 
relation of /i„ oc Aff 6 . 



Figure 6. Stellar radius versus stellar mass for the dEs in the 
range Ah < 3 X 1O 7 M . Best fit is i?* oc Af° 3 . Shown for com- 
parison is R„ oc M2' 2 (see the text). 



rest of the LL family (Fig. This is because the gas that 
remains in dl haloes (or falls back in later on) allows star 
formation to continue after reionization, bringing M* close 
to its value predicted by the energy requirement r\ oc V 2 . 

If radiative feedback is indeed important in the forma- 
tion of dEs, we should try to understand why the scaling 
relations shown by the dEs for the metallicity and surface 
brightness versus stellar mass seem to be natural extrapo- 
lations of those shown by the dls, which we associated with 
supernova feedback. For example, if all dEs have similar 
haloes and therefore similar initial gas masses, and if the 
metals were assumed to be uniformly distributed through- 
out the gas as assumed for larger LLs, then we might have 
expected a steeper dependence of Z oc M*/M oc M* within 
the dE range (rather than the global Z oc M® A which ap- 
plies throughout the LL range). However, since some of the 
gas is expected to photo-evaporate or be kept away from the 
halo even before it cools and falls into the halo center, we 
expect the metals to enrich a smaller fraction of the initial 
gas in fainter dEs, which should lead to a weaker dependence 
of Z on M*, perhaps as flat as Z oc M° A . Given the lim- 
ited width of the dE range, we can probably tolerate some 
deviation from the global Z oc Mj' 4 there. 

As for the surface brightness in dEs, the similarity in V 
between all the dE haloes would lead to the predictions R, ~ 
\R oc \a 1/2 and therefore /i, oc A/*/(Ai?) 2 oc A~ 2 a _1 M* 
(compared to the global /u» oc Mj' 6 ). If the sequence of A/* 
in dEs indeed represents variations in formation time, then 
the combination of the a and A factors should be responsible 
for the flattening of the M* dependence of y,* to the observed 
jU» oc M° . While a is expected to be smaller for larger M» , 
the baryonic spin parameter is expected to be smaller for 
dEs of smaller M», those that formed later and closer to 
Zi on . This is because in those only the gas from the inner halo 
managed to form stars before the reionization time, and this 
inner gas is naturally expected to be of lower than average 



8 DISCUSSION 

We identify four basic characteristic scales in the theory of 
galaxy formation, each originating from a different physical 
process, and each having a different imprint on the galaxy 
population, as follows: 

(i) The upper limit for bright galaxies separating them 
from clusters of galaxies, at M* ~ 10 12 Mq, is where ra- 
diative cooling occurs on a dynamical time scale (Rees & 
Ostriker 1977; Silk 1977; White & Rees 1978 in the context 
of dark haloes). The cooling curve in the relevant temper- 
ature and density range makes this bound roughly coincide 
with an upper bound to the total mass. 

(ii) Supernova feedback becomes effective in heating and 
removing gas in haloes of V < lOOkms -1 , as predicted by 
Dekel & Silk (1986). We argue that this scale marks the 
transition between HH and LL/dwarf galaxies, as seen in 
the SDSS data near M* ~ 3 x 10 10 M Q . 

(iii) Radiative feedback after Z[ on heats the gas to T ~ 
10 4 K, which causes efficient evaporation from haloes of 
V < 20kms _1 (Shaviv & Dekel 2003) and prevents fur- 
ther infall into haloes of V < 30kms _1 (see a review by 
Barkana & Loeb 2001). We propose the possibility that this 
scale, corresponding to A/* ~ 3 x 1O 7 M0, marks the tran- 
sition between the gas-poor dE galaxies and the gas-rich dl 
galaxies. 

(iv) A sharp lower bound for haloes that can form galax- 
ies, at V ~ 10kms _1 , arises from the sharp drop in the cool- 
ing rate below T ~ 10 4 K, where, in the absence of metals, 
it relies on molecular Hydrogen (see Barkana & Loeb 2001, 
Fig. 12). The H2 molecules are dissociated by the weak UV 
flux from the first stars or AGNs long before Zj on (Haiman, 
Rees & Loeb 1996; Haiman, Abel & Rees 2000), allowing no 
gas cooling in these small dark haloes (see DS, Figs. 5,6). We 
find a hint for this lower bound to V in the dwarf spheroidals 
of the Local Group, Fig. [3] 

8.1 HH galaxies 

Therefore, the cooling upper limit and the supernova scale 
limit the stellar masses of bright galaxies to the range 
3 x 10 10 < Ah < 1O 12 M . A significant fraction of the gas 
is assumed to have turned into stars in these galaxies, such 
that 77 is not significantly correlated with the halo proper- 
ties. Then the tight TF relation, the high surface brightness 
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and metallicity, and the weak correlation between the last 
two and stellar mass, all follow naturally from the simplest 
possible assumptions as described in JHl(see also, e.g., Blu- 
menthal el al. 1984). In summary, the assumptions made 
are: 

(i) The halo is in virial equilibrium after spherical collapse 
from a cosmological background. 

(ii) The epoch of galaxy formation is only weakly corre- 
lated with halo mass, consistent with the ACDM cosmology 
where the power index of density fluctuations is n — 2 in 
the range corresponding to HH galaxies. 

(iii) The stellar mass is proportional to the total mass, 
M, oc M, such that rj, the ratio of stellar to initial-gas mass, 
is uncorrelated with the halo mass. This assumption distin- 
guishes the HHs from the LLs. 

(iv) The size of the stellar system is related to the halo 
virial radius by conservation of angular momentum, namely 
via the spin parameter A, which is uncorrelated with the 
halo mass. 

The idealized picture is disk formation via gas contraction in 
pre-formed haloes, unperturbed by recent strong galaxy in- 
teractions. This should be especially valid for the fragile LL 
disks discussed later, where it is supported by a weak spatial 
correlation observed between LSBs and other galaxies, espe- 
cially below a pair separation of 2/t -1 Mpc (Mo, McGaugh 
& Bothun 1994). 

It is not hard to understand why the galaxies at the 
bright end, mostly large ellipticals, are dominated by a 
spheroidal stellar component with little gas and a low cur- 
rent star formation rate (SFR). The high density of the 
cooled gas in the early progenitors of these haloes allows 
the formation of molecules which provide efficient cooling 
even after the gas has cooled to below 10 4 K. This explains 
the high SFR in discs early on. Mergers of discs lead to 
bulges and elliptical galaxies, which therefore tend to be 
those galaxies that dominate the high M, end. The decrease 
in number of objects due to mergers may partly explain the 
low scatter in radius and surface brightness at a given M* in 
the HH regime, as indicated by K03. The mergers provide an 
additional trigger for a high early SFR. The associated high 
gas consumption in these early epochs results in gas-poor 
systems with low SFR today. Since the supernova feedback 
energy is weak compared to the depth of the potential wells 
in HHs, it has a negligible effect on the SFR. 

A potential caveat in the picture that assumes no gas 
loss from HHs may arise from the indications that the bary- 
onic fraction in these systems may in fact be lower than 
the universal fraction by a factor of 2 or more (e.g., Klypin, 
Zhao & Somerville 2002, based on semi-analytic modeling of 
the Milky Way within the ACDM scenario, and references 
therein). Such gas loss may be expected if the HHs result 
from a hierarchical merger process, where the gas is lost 
at early stages from the small building blocks by supernova 
and radiative feedback. Another possibility is that supernova 
feedback is actually stronger than assumed, either due to mi- 
croscopic effects such as porosity in a multiphase ISM, or due 
to hypernova from very massive stars (Silk 2003). Alterna- 
tively, there might be an even stronger feedback mechanism 
at work in big galaxies and in clusters. Hints from SDSS for 
a correlation between HHs and AGN activity (Kauffmann 
el al. 2003, in preparation), together with the established 



presence of massive black holes in early-type galaxies and 
the known energetic radio jets associated with AGNs, may 
provide a clue for the required energetic feedback process. 

8.2 LL and dwarf galaxies 

Most of the galaxies and most of the mass belong to the 
LL and dwarf family below the transition scale: Af* < 
3 x 10 10 Mq. Their halo velocities are below the critical su- 
pernova scale of Vsn ~ lOOkms -1 and they are therefore 
subject to supernova feedback effects which can determine 
their characteristic scaling relations, as argued based on the 
simplest possible model in J^j The energy fed into the gas 
leads to a lower stellar mass fraction M*/M and therefore 
lower surface brightness and metallicity in haloes of lower 
V. Some of the gas may be blown out and some may be 
retained or may fall in at a later time. This gas is kept hot 
and possibly turbulent such that the SFR is regulated by 
supernova feedback, as well as by radiative feedback at the 
lower part of the dwarf sequence. Note that for the scal- 
ing relations to be valid in the LL regime the gas does not 
have to be blown away — it should just be prevented from 
forming stars too efficiently. Our feedback model predicts 
M t /Mg oc V 2 , where M g is the mass of the gas affected 
by feedback and prevented from forming stars. In LLs and 
dwarf irregulars, a significant fraction of this gas must have 
been retained in the galaxy rather than been blown away. In 
this case, the model indeed predicts an increasing gas-to-star 
ratio for decreasing halo mass, as observed. 

Our key assumption for supernova feedback is that 
-Esn oc M». It is crudely justified also in the presence of 
significant radiative cooling, when the gas is at T ~ 10 5 K, 
based on the analysis of supernova remnants by DS. The sec- 
ond assumption is the straightforward energy requirement 
for affecting most of the original gas, -Esn oc M g V 2 . To- 
gether they yield that the effectiveness of feedback varies 
along the LL sequence as 

M t /M oc V 2 . (33) 

The scaling relations for LLs are then obtained using the 
same standard assumptions as used for HHs, namely virial 
equilibrium after spherical collapse and angular-momentum 
conservation, noting that the correlation of formation time 
with halo mass is even weaker for dwarfs where n — * —3. 
Our basic energy condition is clearly based on a simplistic 
model for feedback, which was expected to provide rough 
estimates at best. The fact that this model recovers so well 
the observed scaling relations is partly a matter of lucky co- 
incidences and it should not be taken too literally. However, 
our main moral from the remarkable success of the crude 
model is that supernova feedback can be the primary phys- 
ical process determining the fundamental line of LL/dwarf 
galaxies. One may in fact reverse the logic and infer the 
feedback energy relation, Esn oc M» oc M g V 2 , from the ob- 
served scaling relations, via the other standard assumptions 
of virial equilibrium and spherical collapse in ACDM cosmol- 
ogy. Our toy analysis therefore provides the basis and the 
motivation for detailed future studies of the supernova feed- 
back effects, using more sophisticated modeling and simula- 
tions. The inferred energy relation should serve as a useful 
constraint that must be obeyed by these models. It may be 
a non-trivial challenge for such realistic models to achieve a 



Feedback in dwarf/LSB galaxies 13 



match with observations as good as the match achieved by 
the naive toy model. 

The actual supernova feedback process is likely to be 
much more complex than assumed in our toy model. For ex- 
ample, suprenovae exploding in a disk would affect the disk 
gas and the halo gas in different ways and in an aspheri- 
cal configuration, with fountains punching out the ISM and 
the IGM in a nonuniform and possibly porous manner (e.g., 
MacLaw & Ferrara 1999; Scannapieco, Thacker & Davis 
2001; Mori, Ferrara & Madau 2002; Scannapieco, Ferrara & 
Madau 2002; Silk 2003). This would affect the way the ISM 
and IGM are enriched with metals (e.g., Madau, Ferrara & 
Rees 2001; Thacker, Scannapieco, & Davis 2002). Another 
complication is that the supernova energy can be transferred 
to the gas in either bulk kinetic energy or thermal energy, 
but it can also be kept in reservoirs of other forms such 
as turbulence, which would amplify the feedback effects on 
the gas (e.g., Efstathiou 2000; Thacker & Couchman 2000, 
2001; Springel 2000; Springel & Hernquist 2003a,b). These 
preliminary studies would hopefully guide us to a reliable, 
realistic, detailed feedback model, but our results indicate 
that the global energy balance should somehow mimic the 
predictions of the early naive toy models (Larson 1974; DS). 

If supernova feedback is the primary parameter varying 
along the fundamental line, the other factors affecting the 
galaxy properties can be regarded as secondary parameters, 
responsible for the scatter about the mean scaling relations. 
These may include the following: 

(i) The spin parameter A and the internal angular- 
momentum distribution for the dark matter and in partic- 
ular for the baryons (e.g., Dalcanton, Spergel & Summers 
1997). 

(ii) The halo density profile, e.g., as parameterized by its 
concentration parameter (Navarro, Frenk & White 1997; see 
Bullock et al. 2001a). 

(iii) The epoch when the halo collapsed relative to the 
typical collapse time of haloes of a similar mass and its mass 
accretion/merger history. 

(iv) The detailed efficiency of gas cooling and gas removal 
by feedback. 

(v) The star formation efficiency, e.g., as a function of gas 
surface density. 

In general, the surface brightness at a given M* is expected 
to be below the mean relation for younger galaxies that had 
more recent major mergers and therefore typically have high 
spin and low concentration (e.g. Wechsler et al. 2002). These 
qualitative predictions are consistent with the findings from 
the SDSS data as well as the Local-Group dwarfs. 

As mentioned in the distribution of spin parame- 
ter, in particular, is a natural cause for variation in stel- 
lar radius, and therefore surface brightness of disk galaxies 
(see modeling by Fall & Efstathiou 1980; Dalcanton, Spergel 
& Summers 1997; Mo, Mao & White 1998; van den Bosch 
2000, 2001; Avila-Reese & Firmani 2000, 2003; Firmani C. 
& Avila-Reese 2000). All the scatter observed in R, and y,* 
for SDSS galaxies (K03) can be accounted for by scatter in 
halo spin parameter (e.g., Bullock et al. 2001b). However, 
the lack of correlation between A and halo mass implies that 
the spin parameter cannot be the primary factor driving the 
systematic variation of fi» with halo mass for LL/dwarfs. 

It is worth noting though that the key relation lead- 



ing to the observed scaling relations, M*/M oc V 2 , might 
in principle arise from processes other than feedback. For 
example, M*/M may increase with halo mass because of 
a systematic increase in star formation efficiency due to 
the higher gas surface density (e.g., Galaz et al. 2002). Is 
this enough for explaining the observed correlation without 
a contribution from feedback effects, namely such that the 
cold-gas mass is oc M? Consider a star formation rate which 
depends on gas surface density as in the Kennicutt-Schmidt 
law, /I* oc fig^s, namely M*/M oc fMg^ 1 . From the virial re- 
lation w6 hciv6 from gq. £Lnd 6C[. Mgas oc ■oc 
X~ 2 V 3t with «S r < 1/3 (for -2 > n > -3). We therefore 
obtain M,/M oc ^ 3r (iv-i) _ If 

we assume M* oc M*, as in- 
dicated for the SDSS LLs and the Local-Group dls, we get 
the required scaling relation M*/M oc V 2 for 3t(N — 1) = 2. 
With t ^ 1/3 we obtain N 3. This is a stronger de- 
pendence of SFR on surface gas density than measured for 
star-forming galaxies, N = 1.4 ± 0.15 (Kennicutt 1998). It 
indicates that star-formation efficiency is not a natural pri- 
mary driver for the LL scaling relations, though it may have 
an important role. 

8.3 Related work on the role of feedback 

Our model for the additional role of radiative feedback in 
distinguishing dwarf spheroidals from dwarf irregulars can 
be regarded as a qualitative speculation, to be investigated 
in more detail in future work. 

The trend of increasing total-mass-to-light ratio with 
decreasing mass, which naturally results from supernova 
feedback and reproduces the fundamental line, may not 
be enough for fully resolving the mystery of missing dwarf 
galaxies (Klypin et al. 1999; Moore et al. 1999; Springel et 
al. 2001). The discrepancy is not only between the faint-end 
slopes of the predicted mass function and the observed lu- 
minosity function, but it also involves the distribution of in- 
ternal velocities. The halo masses inferred straightforwardly 
from the observed velocities in the LG dwarf galaxies, many 
of which are on the order of ~ lOkrns -1 , are too small com- 
pared to the CDM model predictions. One possibility is that 
the relevant halo velocities are severely underestimated be- 
cause the sampling by stars in dwarfs is biased towards the 
very inner halo regions, where the rotation and dispersion 
velocities may be significantly lower than the maximum or 
virial velocity relevant for mass estimation (Stoehr et al. 
2002). This possibility is unlikely to provide the full an- 
swer because the velocities measured from HI gas, which 
typically samples more extended radii by a factor of 2-3, 
are still of similarly low, ~ lOkrns -1 (Blitz 2003, private 
communication). The discrepancy between the observed ve- 
locity function and that predicted by the CDM model thus 
seems to indicate the presence of some barren haloes, which 
are completely dark and show no trace of luminous stars in 
them (e.g., Kochanek 2001). Attempts have been made to 
explain the barren haloes by the radiative feedback effects 
discussed in H7.2I which indeed are expected to "squelch" 
the formation of stars in haloes that form after cosmological 
reionization at z ~ 7 (e.g., Bullock, Kravtsov & Weinberg 
2000; Somerville 2002; Tully et al. 2002). Such barren dark 
haloes may alternatively be explained by the destructive 
effect of energetic outflows from one galaxy on neighbor- 
ing forming protogalaxies via ram pressure brushing aside 
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the tenuously held gas (e.g., Scannapieco, Ferrara & Broad- 
hurst 2000; Scannapieco & Broadhurst 2001; Scannapieco, 
Thacker & Davis 2001; Thacker, Scannapieco & Davis 2002). 
These processes are yet to be studied with more realistic 
simulations, in an attempt to find out whether they can in- 
deed explain the absence of luminous components in massive 
enough haloes. 

A scale similar to the supernova scale, originating in a 
different way from the features of the cooling curve, is asso- 
ciated with another transition between different behaviors, 
as discussed in Birnboim & Dekel (2003). They found, using 
analytic arguments supported by simulations with a spher- 
ically symmetric Lagrangian hydrodynamical code, that in 
haloes less massive than ~ 3 x 10 1 Mq the gas falling into 
the halo does not cross a virial shock until it hits the "disc" 
itself. The "standard" virial shock develops only in more 
massive haloes, hosting large galaxies and clusters, where 
the shock quickly expands to near the virial radius. Then, 
as commonly assumed, infalling gas is heated behind the 
shock to the halo virial temperature and is kept pressure 
supported in the halo until it cools radiatively and slowly 
contracts into the disc. In less massive haloes, where the 
virial temperature is below a few x 10 K, the shock that 
tries to develop loses energy very efficiently via radiation 
that is dominated by He recombination and Oxygen lines. 
This prevents the shock from ever expanding into the halo. 
A possible implication of this result is that early star forma- 
tion becomes more efficient in haloes of M < 3 x 10 1 Mq, 
in shocks produced by the cold infalling gas when it hits the 
cold gaseous disc, giving rise to the burst which heats much 
of the remaining gas and produces an LSB galaxy. Further 
infalling gas may prevent blowout and keep the hot gas in 
the galaxy, giving rise to gas-rich dwarf irregulars. 

The low baryonic fraction observed in V ~ 60kms~ 1 
LSB's (by BBS), together with the prediction of r\ oc V 2 
for the feedback effect, implies Vsn ~ 80kms _1 , which is 
quite consistent with the DS estimate of Vsn ~ 100 km s^ 1 
and with the observed transition scale at a stellar mass 
M» £ 10 10 M Q . Furthermore, as argued by Mailer & Dekel 
(2002), feedback can help solve the apparent angular- 
momentum problem within the CDM scenario, where the 
baryons in cosmological simulations seem to lose most of 
their angular momentum and fail to form large discs as 
observed (Navarro & Steinmetz 2000). Mailer, Dekel & 
Somerville (2001) and Mailer & Dekel (2002) modeled the 
properties of the LSB galaxies observed by BBS based on 
spin buildup from the orbital angular momenta along the 
halo merger history, combined with gas blowout from small 
merging satellites and the associated baryonic spin increase. 
They found that a value of Vsn ~ 90kms _1 can indeed 
explain the higher spin observed in the LSB galaxies. 

We note that gas blowout in small haloes may also help 
resolve a third problem of the CDM scenario, the cusp/core 
problem of halo density profiles, where simulated haloes 
show a steep inner cusp while observations indicate that 
at least some galaxies have flat-density cores. While feed- 
back cannot significantly affect the dark-matter distribution 
in big galaxies, an impulsive blowout may reduce the core 
densities in dwarf satellites by a factor of a few (Gnedin & 
Zhao 2002). When these puffed- up satellites merge to build 
up bigger haloes, they get tidally disrupted before they man- 
age to penetrate the inner regions and turn the cores into 



cusps. In this indirect way the feedback can help the sur- 
vival of cores even in relatively big galaxies (Dekel, Devor & 
Hetzroni 2003; Dekel et al. 2003). However, working out a 
feedback model within the CDM scenario that will explain 
the possible existence of cores in giant galaxies and clusters 
of galaxies could be challenging; it will require a feedback 
mechanism more energetic than simple supernova feedback, 
perhaps by hypernovae from massive stars or by radio jets 
from AGNs. 

We conclude that feedback effects seem to be able to 
provide the cure to all three major problems facing galaxy 
formation theory within the CDM scenario. Understanding 
the details of the feedback processes is therefore a major 
goal of galaxy formation studies. 
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